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* 
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Abstract —In this paper, we present distributed controllers for 
sharing primary and secondary frequency control reserves for 
asynchronous AC transmission systems, which are connected 
through a multi-terminal HVDC grid. By using Lyapunov 
arguments, the equilibria of the closed-loop system are shown 
to be globally asymptotically stable. We quantify the static 
errors of the voltages and frequencies, and give upper bounds 
for these errors. It is also shown that the controllers have 
the property of power sharing, i.e., primary and secondary 
frequency control reserves are shared fairly amongst the AC 
systems. The proposed controllers are applied to a high-order 
dynamic model of of a power system consisting of asynchronous 
AC grids connected through a six-terminal HVDC grid. 

I. Introduction 

Transmitting power over long distances with minimal 
losses is one of the greatest challenges in today’s power 
transmission systems. The strong rising share of renewables 
increased the distances between power generation and con¬ 
sumption. This is a driving factor behind the development 
of long-distance power transmission technologies. One such 
example are large-scale off-shore wind farms, which often 
require power to be transmitted in cables over long distances 
to the mainland power grid S). Due to the high resistive 
losses in AC cables, high-voltage direct current (HVDC) 
power transmission is a commonly used technology for power 
transmission in these cases. The higher investment cost of an 
HVDC transmission system compared to an AC transmission 
system, is compensated by the lower resistive losses for 
sufficiently long distances 03 . The break-even point, i.e., 
the point where the total construction and operation costs 
of overhead HVDC and AC lines are equal, is typically 
500-800 km E). However, for cables, the break-even point 
is typically less than 50 km 03. Increased use of HVDC 
technologies for electrical power transmission suggests that 
future HVDC transmission systems are likely to consist of 
multiple terminals connected by several HVDC transmission 
lines 0. Such systems are referred to as Multi-terminal 
HVDC (MTDC) systems in the literature. 

Maintaining an adequate DC voltage is the single most 
important practical control problem for HVDC transmission 
systems. If the DC voltage deviates too far from the nominal 
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operational voltage, equipment could be damaged, resulting 
in loss of power transmission capability and high costs. 

Many existing AC grids are connected through HVDC 
links, which are typically used for bulk power transfer 
between AC areas. The fast operation of the DC converters, 
however, also enables frequency regulation of one of the 
connected AC grids through the HVDC link. One practical 
example of this is the island of Gotland in Sweden, which is 
only connected to the Nordic grid through an HVDC cable 
0. However, since the Nordic grid has orders of magnitudes 
higher inertia than the AC grid of Gotland, the influence of 
the frequency regulation on the Nordic grid is negligible. By 
connecting several AC grids by an MTDC system, primary 
frequency regulation reserves may be shared, which reduces 
the need for frequency regulation reserves in the individual 
AC systems HD- In El, distributed control algorithms have 
been applied to share primary frequency control reserves of 
asynchronous AC transmission systems connected through an 
MTDC system. However, the proposed controller requires a 
slack bus to control the DC voltage, defeating the purpose 
of distributing the primary frequency regulation reserves. 
In ID, distributed controllers for secondary voltage control 
of MTDC systems are proposed, which do not rely on a 
slack bus. The analysis in the aforementioned reference is 
however restricted to the dynamics of the MTDC system, 
thus neglecting any dynamics of connected AC systems. 

In Q and E), decentralized controllers are employed to 
share primary frequency control reserves. In na no stability 
analysis of the closed-loop system is performed, whereas 0 
guarantees stability provided that the connected AC areas 
have identical parameters and the voltage dynamics of the 
HVDC system are neglected. In fl7l . optimal decentralized 
controllers for AC systems connected by HVDC systems are 
derived. In contrast to the aforementioned references, 10 
also considers the dynamics of connected AC systems as 
well as the dynamics of the MTDC system. A distributed 
controller, relying on a communication network, is proposed 
in 0. Stability is guaranteed in the absence of communi¬ 
cation delays. The voltage dynamics of the MTDC system 
are however neglected. Moreover the implementation of the 
controller requires every controller to access measurements 
of the DC voltages of all terminals. 

In this paper we present two distributed controllers for 
combined primary and secondary frequency control of asyn¬ 
chronous AC systems connected through an MTDC system. 


The first controller requires that the communication network 
constitutes a complete graph, while the second controller only 
relies on local information from neighboring AC systems. In 
contrast to existing controllers in the literature, the proposed 
controllers in this paper consider the voltage dynamics of 
the MTDC system, in addition to the frequency dynamics of 
the AC systems. The equilibrium of the closed-loop system 
is shown to be globally asymptotically stable for any set of 
controller parameters. We furthermore bound the asymptotic 
errors of the voltages and frequencies at the equilibrium, 
and show the achievable performance is better than for the 
corresponding decentralized controller studied in jZ|- We also 
show that the frequency control reserves are asymptotically 
shared approximately equally, which is referred to as power 
sharing in the literature. 

The remainder of this paper is organized as follows. In 
Section [IIJ the mathematical notation is defined. In Section 
m the system model and the control objectives are defined. 
In Section [Tv| we recall the decentralized proportional con¬ 
troller for distributing primary frequency control proposed 
in 0. In Section [V] two secondary frequency controllers 
for sharing primary and secondary frequency control reserves 
are presented. In Section [Vi] simulations of the distributed 
controller on a six-terminal MTDC test system are provided, 
showing the effectiveness of the proposed controller. The 


paper ends with concluding remarks in Section VII 


II. Preliminaries 

Let Q be a static, undirected graph. Denote by V = 
{1,..., n} the vertex set of Q, and by £ = {1,..., to} the 
edge set of Q. Let A/} be the set of neighboring vertices to 
i £ V. Denote by B the vertex-edge incidence matrix of Q, 
and let £\,y = BWB 1 be the weighted Laplacian matrix of 
Q, with edge-weights given by the elements of the diagonal 
matrix W. We denote the space of real-valued n x to- valued 
matrices by R nxm . Let CL denote the open left half complex 
plane, and C _ its closure. We denote by c nxm a vector or 
matrix of dimension n x to, whose elements are all equal 
to c. For a symmetric matrix A, A > 0 {A > 0) is used 
to denote that A is positive (semi) definite. I n denotes the 
identity matrix of dimension n. For simplicity, we will often 
drop the notion of time dependence of variables, i.e., x(t) 
will be denoted x. Let ||-|| denote the maximal absolute 
value of the elements of a vector. 


III. Model and problem setup 

We will here give a unified model for an MTDC system 
interconnected with several asynchronous AC systems. We 
consider an MTDC transmission system consisting of n con¬ 
verters, each connecting to an AC system, denoted 1 ,,n. 
The converters are assumed to be connected by an MTDC 
transmission grid. The dynamics of converter i is assumed to 
be given by 

= -±-(Vi-V j ) + 1 ?, (1) 

jeW " 


where V) is the voltage of converter i, C, >0 is its 
capacitance, and J‘ nj is the injected current from an AC grid 
connected to the DC converter. The constant R l7 denotes 
the resistance of the HVDC transmission line connecting the 
converters i and j. The graph corresponding to the HVDC 
line connections is assumed to be connected. The AC system 
is assumed to consist of a single generator which is connected 
to the corresponding DC converter, representing an aggregate 
model of the AC grid. The dynamics of the AC system are 
given by the swing equation El: 

mM = Pf n + P™ m + P™ - Pf, (2) 

where TOj > 0 is its moment of inertia. The constant P J nom 
is the nominal generated power, PJ" is the uncontrolled 
deviation from the nominal generated power, P‘ n] is the 
power injected to the DC system through the convertera, P® en 
is the generated power by the generation control (primary or 
secondary) of generator i, respectively. The control objective 
can now be stated as follows. 


Objective 1. The generation control action should be asymp¬ 
totically distributed fairly amongst the generators, i.e.. 


lim 

t—¥ OO 


Pf en {t) + 


1 

n 


Y. p i 


< e ge " 


Vi = 1,..., n, 


where e gen is a given scalar. Furthermore, the frequencies of 
the AC systems, as well as the converter voltages, should not 
deviate too far from their nominal values, i.e., 

lim | Vi(t) — V[ ef \ < e l Vi = 1,... ,n 

t—> OO 

lim |w,;(t) — u) ref \ < e“ Vi = 1,... ,n, 

i—>■ oo 


where Vf '^ is the reference DC voltage of converter i, u> re f is 
the reference frequency and e v and e“ are given scalars. 

Remark 1. It is in general not possible to have e v = 0, 
since this does not allow for the currents in the HVDC grid 
to change. The bounds e gen and e“ on the other hand, can 
in theory be zero. However, this may be hard to achieve with 
limited communication. 

IV. Decentralized MTDC control 

In this section we summarize previous results on de¬ 
centralized control of MTDC systems. The detailed results 
and proofs can be found in |2). We consider the following 
decentralized frequency droop controllers for the control of 
the AC systems connected through an MTDC network 

= (3) 

where w, is the frequency of the generator, u/ ef is the 
reference frequency (e.g., 50 or 60 Hz), and /f‘ lroop > 0 . 
The local controllers governing the power injections into the 
MTDC network are given by 

pmj = pinj, nom + _ ^ref) + R V (^ref _ (4) 





where Pf’ nom is the nominal injected power, and K'f > 0 
and ATf > 0 are positive controller gains for all i = 1 ,..., n. 
The HVDC converter is assumed to be perfect and instanta¬ 
neous, i.e., injected power on the AC side is immediately 
converted to DC power without losses. Furthermore the 
dynamics of the converter are ignored, implying that the 
converter tracks the output of controller (|4]> perfectly. This 
assumption is reasonable due to the dynamics of the converter 
typically being orders of magnitudes faster than the AC 
dynamics ifTOll . The relation between the injected HVDC 
current and the injected AC power is thus given by 

V t ^=F». (5) 

By assuming V, = V Dom Vi = 1 ,,n, we obtain 

^nomjinj _ pin] 


Furthermore, we assume that the nominal generated power 
equals to the nominal injected power. 

Assumption 1. Pf m = Pf' nom Vi = 1, ..., n. 

Theorem 1. The equilibrium of the decentralized MTDC 
control system with dynamics given by 0. 0. 0. 0, and 
where 0 holds, is globally asymptotically stable /|2]i. 

We will now study the asymptotic voltages and frequencies 
as well as the generated power of the MTDC system. We 
make the following assumption on the controller gains. 


Assumption 2. The controller gains satisfy K'f = 
k u , Kf oop = k dmop , KY = k v Vi = 1,..., n. 

Theorem 2. Given that Assumption [2] holds, then for the 
HVDC and AC systems 0. 0 with generation control 0 
and converter control 0 and where 0 holds, Objective [7] is 
satisfied for e gm = e 8 d e " c , e v = e dec and e“ = e dec 0, where 


k dm °P max, Pf 


dec pdroop _j_ fouj 

kP 


k v ” 

(n — 1 ) + y 


V nom ^ A i(C R ) 


P v - 
°dec 


e dec = 


nk dm °Pk v 
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Ear 

»=t 


kP max.j | P” 

{kP +k dr °op)V" om A i{C R ) 


E 


nk dm °P 




i=i 

maxi Pf 

k dm °P + k u 


(n-1 ) + t^E 


V nom A i(C R ) j ' 


In the following section we will show that the above upper 
bounds can be tightened by a secondary control layer. 


V. Distributed secondary frequency control 

In this section we study a distributed secondary frequency 
controller to tighten the stationary error bounds of the closed- 
loop system given in Theorem [2] To implement a distributed 
controller, we assume the existence of a communication, 
enabling communication between all AC generators. We will 
consider two distributed secondary controllers; a controller 



(a) (b) 


Figure 1: Illustration of the HVDC grid and the communi¬ 
cation network topologies. The HVDC lines are illustrated 
with solid lines, while comminication lineas are illustrated 
with dashed lines. 0 shows the structure of the complete 
communication structure of 0 , while 0 shows the structure 
of the communication structure of 0 . 


where any pair of generators can communicate directly, and 
a controller where only neighbouring generators can com¬ 
municate directly. These two controllers correspond require 
a complete communication graph and a communication graph 
containing a spanning tree, respectively. The corresponding 
topologies of these controllers are illustrated in Figure [T| The 
first generation controller of the AC systems is given by 


TSV 1 n 

pf „ = _^oop (wi _ wref) _ cop, IA £ , 


7^droop,I / ref \ 

Vi = Ki [Ui - id ) - 777 , 


(7) 


where 7 , Af oop ' \ i = 1 are positive constants. If 

7 = 0 , then 7 ,; becomes a scaled integral state of the 
local frequency deviation (u 7 — a/ ef ). The second generation 
controller of the AC systems is given by 

TfV 

Pf n = -Af oop ( Wi - w ref ) - gyAf oop ' 7 
Vi = ^f 00p ' I (wj - W ref ) - 6 ^ djirji - rtj). (8) 

jeNi 


We assume Cij = c :)l , i.e., the communication graph is 
undirected. The above controller can be interpreted as a 
distributed Pi-controller, with a distributed consensus filter 
acting on the integral states //,. In contrast to the controller 
0, 0 does not rely on a complete communication graph. 
Note that we make no assumption that the communication 
topology resembles the topology of the MTDC system. 


A. Stability analysis 

Combining the voltage dynamics 0 , the frequency dy¬ 
namics 0 and the generation control 0 or 0 . the 
converter controller 0 and the power-current relation¬ 
ship 0 , together with Assumption [T| and defining 7 = 
[ 71 ,..., 7 n] r , Cj = uj — w ref l„ and V = V — V lef , 
where V ref = [Vj ef ,..., Vf f ], we obtain the closed-loop 
dynamics given by 0 or m respectively where M = 
diag(mi _1 ,..., m„ _1 ) is a matrix of inverse generator 
inertia, E = diag([Cf 1 ,..., Cf 1 ]) is a matrix of electrical 
elastances, A " droop - 1 = diag([A'f°° P '\ ..., A'*°° p ’ 1 ]), and C c is 
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the weighted Laplacian matrix of the communication graph 
with edge-weights c,j. The details of the derivations have 
been omitted, but the derivation follows the steps taken 
in Section IV in |2). Clearly the representation ([9]) is not 
minimal with respect to the output y = [Cj t ,V T ] T , since 
the integral states p are redundant. It is easily shown that 
substituting the projection 


r{ = -l 7 r 


( 12 ) 


in & the output dynamics remain unchanged. Even though 
the dynamics |9]) are redundant, they show that the secondary 
control layer can be implemented distributively. Each AC 
generator needs only to compute its local integral state, 
and communicate this to the remaining generators. In order 
to prove stability of the closed-loop system, it is however 
essential to consider the reduced system after applying the 
projection ( fl2l i, which is given by ©■ This system can be 
interpreted as a centralized implementation of (|9j, where the 
single integral state rf is computed centrally with access to 
all frequency measurements. 

Assume that the system matrices of © and © are 
full-rank, which ensures that unique equilibria of (lOi and 
© exist. Denote these equilibria 2o,i = [<jJq\, 
and i 0 2 = [tU(f 2 , Vo 2 ] T ’ respectively. Define x\ = 
[of, V T, fjf] T = (Cj t ’v t , y T } T Ki, V^nl if and x 2 , 
mutatis mutandis. Now: 


X\ = Ax i (13) 

x 2 = Ax 2 (14) 

with the origin as the unique equilibria of both above 
dynamical systems. We are now ready to show the main 
stability result of this section. 

Theorem 3. The equilibria of the systems defined by © 
and © are globally asymptotically stable. 

Proof: First consider the Lyapunov function candidate 

i T/nom 

W{Q, V, if) = -cj T K UJ (K v )- 1 M~ 1 cD + V T EV 

+ \(V '?■ (15) 


Clearly W(u, V,rf) is positive definite and radially un¬ 
bounded. Differentiating © with respect to time along 
trajectories of (JT3J, we obtain 


W(io,V,fj') 

= u) T K u (K v )~ 1 M~ 1 u> + V nom V T EV + fj'fj' 

= ui T ( - K U (K V )~ 1 (K U + K dmop )uj 

+ K U V - -A" dr0 °P’ h n fj') 
n ' 

+ V t (k u ’Cj - ( V Dom C R +K v )V ) 

+ 77' T (il^/C d ™P■ I cu'-7f? , ) 

= -w T ( - K U ’(K V )~ 1 (K U1 + K dmop )Ci 
+ 2 lu t K u V - V T (V aom C R + K v )V - 7(7)') 2 


- u T V T 

— 7(t? / ) 2 - 


K LJ (K V )~ 1 (K U) + AT droop ) -K u 
—K u K v 



UJ 
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=Q i 


Clearly W(at, V, ff) < 0 iff the symmetric matrix Q\ is pos¬ 
itive definite. By applying the Schur complement condition 
for positive definiteness to Q\, we see that Q\ is positive 
definite iff 


K U (K V )~ 1 (K U + K dmop ) - K U (K V )~ 1 K U 
= (K v )~ 1 K dloop > 0. 


Hence Qi is always positive definite. If 7 > 0, 
W(uj,V,fj') < 0, and the origin of ( fT3j ) is thus globally 
asymptotically stable. If however 7 = 0 then W(Q, V, r}') < 
0 and the set where W(Ci,V,fj') is non-decreasing is given 
by 

G={(Q,V,rf)\W(Q,V,rf) = Q} 

= {(u,V,fj')\rf = k}, 


for any k £ K. Clearly the largest invariant set in G is the 
origin. Thus, by LaSalle’s theorem for global stability, the 
origin of ([T3]l is globally asymptotically stable also for 7 = 0. 



































Consider now the following Lyapunov function candidate 


= e dec and e “ = e dist’ where 




i T/nom 

-uj t K u (K v )- 1 M~ 1 uj + -^-V T CV 

+ lv T fj, (16) 


where C = diag([C'i,..., C n }). Clearly W (w, V, fj) is posi¬ 
tive definite and radially unbounded. Differentiating ( p~ 6 l > with 
respect to time along trajectories of (fl4|. we obtain 


W(u,V,fj) 

= Q t K u (K v )- 1 M- 1 u> + V nom V T EV + fj T rj 
= cd T (- K u (K v )~ 1 (K u + A' droop )w 
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since Q i is positive definite. The set where W(oj,V ,rj) is 
non-decreasing is given by 


G = {(u,V,r])\W(u;,V,rj) = 0} 
= {(u,V,f })\ry = kl n }, 


for any k £ K. Clearly the largest invariant set in G is the 
origin. Thus, by LaSalle’s theorem for global stability, the 
origin of 0 is globally asymptotically stable. ■ 


B. Equilibrium analysis 


In this section we study the properties of the equilibria of 
0 and 0 , respectively. We show that by employing the 
aforementioned secondary frequency control schemes, it is 
possible to tighten the error bounds in Theorem[2] Analogous 
to Section IV we assume uniform controller gains. 


Theorem 4. Assume that A ssumptions [7] and [2] hold, mutatis 
mutandis. Consider the HVDC and AC systems © 0 with 
generation control 0 and where the relation 0 holds. 
Objective [7] is satisfied for the secondary controller 0 when 
7 —> 0 + and for the secondary controller 0 when S —> +oo. 
In both cases, Objective [ 7 ] is satisfied for e dlst = e^" t , 
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V nom E X^Cr) I ' 


Furthermore 1 Jyh = l^V = 0, i.e., the average frequency 
and voltage deviations are zero. 


Proof: Before studying the equilibria of ( fTO} and 0 , 
we will show that under the assumptions that 7 —> 0 + , and 
S —> Too, the first 2n rows of the equilibria of ( |10[ i and 0 
are identical. Consider the last row of the equilibrium of ( fTO} , 
which as 7 —> 0 implies l^Co = 0. Thus, premultiplying the 
(n + l)th to 2 nth rows of the equilibrium of ( fTO} with l^C 
yields 1„V = 0. Finally, premultiplying the first n rows of 
( fTO! with l^M -1 yields 


nk v k dm ° v ' 1 


rf ' = -1 T n P n 


(17) 


Now consider the equilibrium of 0 - Following similar steps 
as in the manipulation of ( p~ 0 ] >, we obtain that = 0 
and l^V = 0. Thus premultiplying the first n rows of the 
equilibrium of (fTT]) with we obtain 


kVkfr oop, 1 


-l 1 n = -l 1 P n 


(18) 


We write 77 as a linear combination 77 = J27=i a i v i’ where 
is the (normed) ith eigenvector of C c . Inserting the eigen 
decomposition of 77 in ( fl 8 | yields 

o = _l W 

7 . V ^.droop, I ’ 

since v? = !—1 ,,. In order to determine a? for i > 2, we 

consider the last n rows of the equilibrium of ( fTTj ). Again, 
using the eigen decomposition of 77 we obtain 


k dr °° P' ^ = SC c = E 

i—1 i—2 


;multiplying the above equation with ( 77 ° ) T we obtain 

7 .droop, I ( V 0)T£ 
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Clearly a,j —> 0 as S —> 00 for j = 2,..., n, if Cj is bounded. 
But Cj must be bounded since the system matrix of (111 
is full rank, implying that the steady-state solution to (111 
is bounded. Thus, at steady-state we have 77 , = 77 * \/i = 
1 Inserting this in (p~ 8 ]> yields 


nk v k imop ’ 1 
AT 


77* = -1 T n P m . 


(19) 


Comparing ( fl~9l ) with ( | 1 1\ , it is clear that the first 2n rows of 
the equilibria of (jTO]) and 0 are identical, and thus define 


























the same solutions. We thus proceed only considering the 
equilibrium of ([T()|, whose last row implies 
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After some algebra, the following expression is obtained 
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Again, write C2> = ]TA =1 a i v i> where vj is the eigenvector of 
A 2 with eigenvalue Af. Let Af denote the eigenvalue of A 2 
with eigenvector -f=l n . By direct computation 


Af = 


fcdroop fc ^( fc droop,I) 


■If 2 
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|/nom ynomfeoj ry ’ 


and clearly Af —> 00 as 7 —> 0 + , implying that af —> 0 
as 7 —► 0+. For the remaining eigenvalues, i.e., i > 2, we 
obtain after some calculations 


k droop + k u 


Af > 


Write V = J 2 i=i a i v h where uf is an eigenvector of Ai 
with the corresponding eigenvalue Af. It is easily verified 
that l rt is an eigenvalue of Ai, which we denote v{. 
Since A\ is symmetric, its eigenvectors can be chosen to 
form an orthonormal basis. By premultiplying ( |22| > with 
(uj) T , and keeping in mind that Aivj = Afttj, we obtain 
dj = (vj) T P m / A]. By direct computation we obtain 

uV I 7. V (/.droop, I\2 \ 

a: = w ((^)+ pyp ). 

by which we conclude that Af —> 00 as 7 -> 0 + . Thus 
ai —> 0 as 7 —> 0 + . For i > 2 we obtain after some calcu¬ 
lations Af > yi (kP + A Thus, we obtain the 

following bound on V: 
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By premultiplying ([23]) with (af), we obtain 
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Premultiplying the first n rows of ( |2Tj ) with pr(CR+y^I n ), 


where we have used the fact that the eigenvectors of Ai 
are also eigenvectors of Cr. Finally, we consider the power 
output of the generation control. Note that lim 7 _ ) ,o ^f = 





















































VI. Simulations 
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Noting that i> 2 = -^=, we have that l^V 2 = 0 for * > 2. By 
letting 7 —>■ 0, the above equation simplifies to 
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Using the previously derived lower bound on A 2 for i > 2, 
we obtain the following bound on the generated power 
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Remark 2. The upper bounds on the AC frequency and the 
DC voltage errors, i.e., e u and e v in Theorem [7] are lower 
than the corresponding bounds in Theorem f2Jgiven the same 
netw’ork and controller parameters. In particular 
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In this section, simulations are conducted on a test system 
to validate the performance of the proposed controllers. The 
simulation was performed in Matlab, using a dynamic phasor 
approach based on go The test system is illustrated in 
Figure [2] The parameters of the MTDC grid are given in 



Figure 2: MTDC test system, consisting of a 6-terminal 
MTDC grid. Each terminal is connected to an IEEE 14 bus 
AC grid, sketched as octagons. The dashed lines illustrate the 
topology of the communication grid of the controller (|8|. 

[I] and are chosen uniformly for all VSC stations. Note that 
we in the simulation also consider the inductances L t j and 
capacitances C,; 7 of the HVDC lines. The capacitances of the 
terminals are assumed to be given by Ci = 0.375 x 10~ 3 
p.u. The AC grid parameters were obtained from oa. 
The generators are modeled as a 6 th order machine model 
controlled by an automatic voltage controller and a governor 
m. The loads in the grid are assumed to be equipped with 
an ideal power controller. The controllers (f3|, (f7| and 


Table I: HVDC grid line parameters 


(*, j) 

Rij [p.u.] 

Lij [10 3 p.u.] 

Ci) [p.u.] 

(1,2), (1,3), (2,4), (3,4) 

0.0586 

0.2560 

0.0085 

(2,3) 

0.0878 

0.3840 

0.0127 

(2.5), (4,5) 

0.0732 

0.3200 

0.0106 

(2,6), (3,5), (5,6) 

0.1464 

0.6400 

0.0212 


E p <" ■ 

i—1 

However, as the bounds are conservative, no conclusion 
about the actual control errors can be drawn. 


were applied to the aforementioned test grid. At time t = 1 
the output of one generator in area 1 was reduced by 0.2 p.u., 
simulating a fault. The communication network of controller 
is illustrated by the dashed lines in Figure [2] Figure [3] 
shows the frequency response of all AC grids for the three 
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Table II: Controller Parameter 

K“ K v K imo p i<r droop ’ 1 -y <5 

lit l i 

9000 110 8 10 0 5 


controllers considered. Figure [4] shows the DC voltages of the 
terminals. Figure [5] shows the total change in the generated 
power within each AC area. It can be noted that immediately 
after the fault, the frequency at the corresponding AC area 
drops. The frequency drop is followed by a voltage drop 
in all terminals, and a frequency drop at all AC areas. The 
frequencies and voltages converge to new stationary values 
after approximately 30 s. We note that the asymptotic error of 
the frequencies and voltages are significantly smaller when 
the controllers 0 and 0 are employed, than when the 
decentralized droop controller (|3]i is employed. The generated 
power is shared fairly between the AC areas. 



Figure 3: Average frequencies in the AC areas for the 
controllers 0.0 and <(8j, respectively. 



Figure 4: DC terminal voltages for the controllers (|3j, ((TJ 
and ([8]), respectively. 


VII. Discussion and Conclusions 

In this paper we have studied controllers for sharing 
primary and secondary frequency control reserves in asyn- 



Figure 5: Total generated power in the AC areas for the 
controllers 0. 0 and 0, respectively. 

chronous AC systems connected through an MTDC system. 
We have reviewed a decentralized droop controller, and 
later expanded this to two distributed secondary frequency 
controllers. The distributed controllers use both local and 
neighboring frequency measurements of the AC grids, as 
well as the local DC voltage measurements. The resulting 
equilibria of the closed-loop system was shown to be globally 
asymptotically stable by using Lyapunov arguments. We also 
showed bounds for the asymptotic deviations of the DC volt¬ 
ages and the AC frequencies from their reference values. The 
obtained bounds are lower than the corresponding bounds 
when using only decentralized droop control. Furthermore 
the generated power from the primary frequency control 
is approximately shared fairly between the AC areas, and 
the error from fair power sharing is bounded. We have 
furthermore demonstrated our results on a 6 terminal MTDC 
system with connected AC systems. Future work will focus 
on eliminating the static errors in the frequencies. 
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